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ABSTRACT 

This paper examines the control of flexible structures for those systems with actuators that are 
modeled by second order dynamics. Two modeling approaches are investigated. First a stability 
and performance analysis is performed using a low order finite dimensional model of the structure. 
Secondly, a continuum model of the flexible structure to be controlled, coupled with lumped 
parameter second order dynamic models of the actuators performing the control is used. This 
model is appropriate in the modeling of the control of a flexible panel by proof-mass actuators as 
well as other beam, plate and shell like structural members. The model is verified with 
experimental measurements. 

INTRODUCTION 

This paper presents a study of active vibration suppression in flexible structure using actuators 
with second order dynamics. First, a low order model of the structure is used to investigate 
stability properties. It is then shown that the common practice of maintaining stability by using 
relative velocity feedback (i. e. the difference between the structure's velocity and the actuator mass 
velocity) does not necessarily lead to the best closed loop performance. 

In addition to a finite dimensional analysis of the effects of actuator dynamics on active vibration 
suppression, an infinite dimensional model is suggested. This model proposes adapting an 
approach presented by previous authors on combined dynamical systems. This approach is 
adapted here to include a control law acting between a (finite dimensional) second order actuator 
and a structure defined by a second order in time partial differential equation. The specific case of 
a beam modeled as an Euler-Bemoulli equation with both internal (Kelvin-Voight) and external 
(air) damping included and actuator is presented. Modal equations are presented in terms of the 
Green's function without actuator dynamics. The case of velocity feedback of the structure at the 
point of attachment on the beam is derived. The open loop equation with point dynamics are then 
verified experimentally. A short list of acronyms used in this paper follows in the appendix. 

SDOF/PMA 

The linear proof mass actuator (PMA) is a solenoid-like device^). Current flowing in a coil of 
wire attached to the frame of the PMA, combined with permanent magnets fixed to the proof mass 
produces an electromagnetic force that accelerates the proof mass. Furthermore, this 
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electromagnetic force produces a reaction force on the PMA frame that can be used for control law 
actuation. By regulating the current supplied to the coil one can control the force produced by the 

actuator. A linear variable differential transformer (LVDT) is used to measure q, the proof mass's 

position relative to the PMA frame. This signal can be differentiated to give fj, the proof mass's 
velocity relative to the frame. An accelerometer attached to the PMA gives the structural 
acceleration, which is integrated to give its velocity, x s . These three measurements provide three 
local feedback paths that can be used for output feedback control. The combination of the proof 
mass and these feedback paths can be modeled as a second order system as shown in fig. 1. 

In fig. 1 the structure to be controlled is represented by a damped single degree of freedom system, 
(M s , K s , C s ), the actuator is also modelled by a single degree of freedom system with an 
accompanying force generator. The spring stiffness, k a , represents the servo stiffness of the 
actuator. The damper, c a , represents the friction inside the actuator. A force generator, f g , is used 
to show the use of either velocity feedback path. The dead mass, md, associated with the actuator 
is fixed to M s . The equations of motion for this system are 


r^?Mo%°l*4 K o °lx4 0 } 

L m p m p J L u c a J u % J L r gJ 


X=[ x s T) ]T, 


Mc=M s + m<i 


( 1 ) 


In order for any control formulation to be useful it must be stable. Therefore, the effects on 
stability of each of the feedback paths must be examined. 

The relative position feedback gain, k a , must be positive^) for the stiffness matrix to be positive 
definite. If k a were zero the actuator would have an uncontrolled rigid body mode. Therefore, a 
first requirement for stability is that k a be positive. 


With relative velocity feedback, i.e. fg = - Cg q, the damping matrix becomes 

( 2 ) 


Which is symmetric and positive definite for all values of Cg >0, making this a stable control law. 
Yet, the term c a +Cg in the damping matrix casts this type ofcontrol into the same formulation as a 
traditional vibration absorber problemf 1 * 10 ). 

If the control force is made proportional to the structure's velocity, i.e. fg=-gx s ,the damping 
matrix becomes 



Which is asymmetric, such that the system 

M X + CX + K X = 0 (4) 


: r c s o ■ 
L 0 Ca+^g . 
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becomes potentially unstable as g increases. Therefore, the question arises: why use fg=-gx s if 
this type of control can become unstable? If one were designing a colocated velocity control law 
that ignored the actuator dynamics, this is exactly what one would do( 9 ). Our purpose here is to 
determine whether using non-colocated velocity feedback will provide performance benefits over 
the use of only colocated, while including the actuator dynamics in the analysis. 

Finally, a combination of both relative velocity feedback and direct velocity feedback is considered 
rather than each path individually. That is f g = - gx s - Cgfj which produces no new stability 
problems, yet provides for the interaction of both feedback paths. In other words, the use of the 
colocated feedback can be used to stabilize the non-colocated feedback. Therefore, the problem 
becomes that of finding an optimal choice for k a , c g , and g. This problem will be formulated as a 
linear quadratic regulator and solved as a parameter optimization. 

The cost function is chosen to be the infinite time integral of the square of the structure's position, 

oo oo 

J f 

q= [ x s ti x s f|]T 

o = I 1 i=j=1 

^1J t 0 otherwise 



The equations of motion are written in state space form 


q=Aq 


A _f 02x2 12x2 1 

A 1 -M-lK -M-lD J 

M=T M r P Sl 

L m p m p J 

n=T Cs 0 1 

1 g Ca+Cgj 

*[*i] 


The cost function can be rewritten as 


J=qT(0)Pq(0) (7) 

where P represents the solution of a Lyapunov equation 

ATP + PA = -Q (8) 

Note, that in comparison to a typical LQR optimal control problem, that is 


J=Jq T Qq + u t Ru dt (9) 

there is no penalty on the control here, i. e. R=0. In fact there is no control, u, in this problem, 
only feedback gains which are treated as parameters. Hence, while this problem resembles an 
optimal control problem it is really a parameter optimization. 
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The cost function J is a function of both the initial conditions and P, where Pisa matrix whose 
elements are a function of the feedback gains. To remove the dependence of J on q(0) the initial 
conditions are treated as a random vector, and the expected value of J is taken as 

E{J} = E{qT(0)Pq(0)} = tr(PZ°) (10) 


Z° = q(0)qT(0) 


( 11 ) 


Z° is a normalized second order moment matrix of the initial conditions. For this example it is 
assumed that the structure and actuator are at rest with the structure given an initial deflection and 

the proof mass at its equilibrium position relative to the structure. Therefore Z° is 


7 o = ( 1 i=j=l 

\ 0 otherwise 

The final expression for the cost function is 

J = Pn (12) 

A necessary condition for optimality is that 


dJf n dPn 
dF F* U dF 


F* 


(13) 


where F represents the set of feedback gains. The matrix P was computed algebraically, wit h the 
aid of MACSYMA, and the required derivatives calculated analytically. With these analytical 
expressions the optimal feedback gains were determined numerically. A special case of this 
formulation is for zero structural velocity feedback, g=0, the equations simplify substantially and 
the results of Juang( j ) are obtained. 


Example 

Consider attaching a PMA to the tip of an undamped cantilevered beam, with the desire to attenuate 
the first bending mode of vibration of this beam. The material properties of the beam are given in 
Table l($. The mass properties of the PMA are given in Table 2. 

The first natural frequency for the beam plus dead mass calculated from a finite element analysis is 
12.66 rad/sec (2.02 Hz). The equivalent bending spring stiffness is calculated to be 158.2 N/m, 
and the structural first modal mass plus actuator dead mass is 0.987 kg. Table 3 gives the optimal 
parameter settings for two conditions: 1) zero structural velocity feedback (i. e., a vibration 
absorber), and 2) with structural velocity feedback and the relative velocity feedback gain held to 
17.5 N-s/m. 


Fig. 2 shows the structural response of both systems to the same initial condition. Note that, the 
settling time, and maximum overshoot of system 2 is superior to that of system 1. Fig. 3 shows 
the response of the actuator mass for this simulation. The major disadvantage of this control law 
can be seen in this figure, which is the actuator of system 2 requires a greater stroke length. 

There are several lessons that can be learned from controlling a single degree of freedom structure 
with a PMA that can be applied to the vibrational control of both multiple degree of freedom and 
distributed parameter structures. Control laws that ignore actuator dynamics may result in closed 
loop instability. The use of only safe or nondestabilizing feedback paths may not yield the best 
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performance. Furthermore, using only relative position and velocity feedback results in a control 
law that is no different than that of a traditional vibration absorber. This type of design tends to 
require low feedback gains, such that the motion of the proof mass is unimpeded. Finally, better 
performance is achieved with structural velocity feedback combined with relative velocity 
feedback. In fact, a high structural feedback gain can only be tolerated in the presence of a high 
relative velocity feedback gain. 

MDOF/MPMA 


Fig. 4 shows two proof mass actuators attached to four degree of freedom structural model. 
Assume that the actuators are not interconnected. That is, any measurements made by PMA one 
are not shared with PMA two and vice versa. In this case we have the same feedback paths as in 
the SDOF/PMA, case just more of them. 

The relative position feedback gain, k p i, must be positive for each actuator i. This requirement 
eliminates the uncontrolled rigid body mode of the actuator. Furthermore, k p j > 0 is necessary so 
that the stiffness matrix for this system is positive definite. Similarly, the combined relative 
velocity gain for each actuator, (c a + c g )i must be positive. The intuitive proof for these 
requirements is that the addition of a damped single degree of freedom to a stable damped multiple 
degree of freedom structure results in a stable system. 

If structural velocity feedback is used, the stability of the system can be extunined using the well 
known Routh-Hurwitz criteria. For high order systems this test becomes difficult to implement. 
An alternative is to examine the system damping matrix. For multiple degree of freedom systems 
that can be described by the system of equations, Eq. (4), stability is guaranteed if the symmetric 
portion of the damping matrix is positive definite, provided that the mass and stiffness matrices are 
symmetric positive definite. The damping matrix C is made asymmetric by the introduction of the 
feedback gains gj in some of the off diagonal terms. Any asymmetric matrix can be written as the 
sum of a symmetric and a skew symmetric matrix. 

C = Csym + G, C S ym = C sym T G = -G"^ 

Therefore, the stability of the MDOF/ MPMA system is guaranteed if the matrix C sym is positive 
definite^). This test can be applied to the system of the previous section. The matrix C S ym for this 
system is 



£ 

2 

c a +Cg 


(14) 


and will be positive definite if the following relations are satisfied 

C s >0 (15) 

4 Cs (c a + c g ) > g 2 (16) 

Determining whether high order matrix with a small number of feedback gains is positive definite 
can become just as tedious as the full Routh-Hurwitz test. Therefore, it is useful to simplify this 
result by examining the definiteness of the principal minors of Csym- The matrix C S y m can be 
written as the sum of a positive definite matrix C c and a sparse matrix of zeros and 2x2 symmetric 
blocks containing gi placed along the diagonal according to the actuator placement. When this 
sparse block diagonal matrix is added to Q to form C S ym the only minors of C c that are affected 
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are the 2x2 blocks added at the position of the actuator coordinates. Therefore, only those minors 
changed by the addition of feedback need be checked. This test leads to the following result for 
each actuator^ 8 ). 


4(C$i + C s j+i)(ca + Cg)i > gi 2 


(17) 


It should be recalled that these results if satisfied ensure stability, but if violated do not imply 
instability. Therefore, we feel that this stability criteria is of a conservative nature. 

DPS/PMA 

In this section we examine a beam modelled as an Euler-Bemoulli beam with both air damping and 
Kelvin- Voight internal strain rate damping. A proof mass actuator is attached to the beam at some 
point, fig. 5. The analysis here follows that of Bergman( 2 ). The equations of motion for this 
system are: 


3^ 34 * 

putt + (Cl + C 2 I— r) Ut + El — u = £ { [Fj(t) - fgi(t) - mdiUtt(hi,t)] 6(x-hj) } 
dx 4 dx 4 j = i 

mpfli + m p iu t t(hi,t) + c p if|i + k p jT|i = f g i(t) 

Fi(t) = -mpiili - m p iu t t(hj,t) + f g i(t) 

fgi(t) “ - Cgj Ut(hj,t) 

r = the number of actuators 


(18) 

(19) 

( 20 ) 
( 21 ) 


Where Cgj is the structural velocity feedback gain. Note that in this case only feedback of the local 
velocity is considered. 

These equations can then be nondimensionalized according to the following definitions 

6(1 - 6i) = I6(x-hi), 

FKt)=§Fi(t), 
fi(0=§fi(t), 

l 3 flCgj 

c g i - 


— u 
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W’ 
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ci 
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s 

II 
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£ 
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The nondimensional equations of motion are, where the overbars have been dropped 
d^ 34 r 

Utx + (£1 + e 2rrr) u t +rrr u = 5^ { [■ Ppi’Hi ■ Ppi u tt(6i»x) - Pdi u Tt(6i>x)] 8(^-5i) j (22) 

3s 4 os 4 i=l 


f|i + uxx(5i,t) + eifji + ctoi^i = - — C g i u x (5i,T) 

Ppi 


(23) 
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The solution to this problem is assumed to be of the form 


u&T) = X a nC0 X n(^) < 24 > 

n=l 

oo 

Tl(t)= X a n( T ) A in X n(5i) 
n=l 


where 


A in — ' 


OCrf 


a 


01 


a n 


4’ 


Ain - Aj n + 


Ctpi 4 

a 0 i 4 - a n 4 


(26) 


and X n (^) are the eigenfunctions of the undamped system. a G i 4 represents the square of the 

actuator's undamped natural frequency, and a n is the eigenvalue corresponding to X n (^). 

For the special case of zero actuator dead mass, the temporal coefficients satisfy 

8n(x) + (ei + £2 a n 4 ) a n (x) + a n 4 a n (x) 


oo 


+ X{c gi O " 

m=l 



-) + e 

m 


«n 4 «m 4 

Oo 8 


- ei 


- E2a n 4 


) 


M-piAinAimX n (^)X m (^) an(x) — 0 


(27) 


The eigenvalues a n 4 are calculated from 


r 


i=l 


«oi 4 «n 4 

a 0 i 4 - a n 4 


G(8 i ,5 j ;a n ) 



0 


(28) 


where G is the Green's function for a clamped - free Euler-Bemoulli beam. These equations are 
then useful for stability analysis, control gain formulation, and experimental verification. 


EXPERIMENTAL RESULTS 

A combined lumped parameter zero gain distributed parameter experiment was performed to verify 
Eq. (28) for the special case of |!di = 0. A clamped-free Euler-Bemoulli beam, whose material 
properties are given in table one, was appended with two identical passive spring-mass systems. 
One was placed in the middle of the beam, x = 1/2, and one at the free end, x = 1. The mass was 
measured to be 49.2 x 10' 3 kg. The spring stiffness was 858 N/m. 

In a test, an accelerometer was placed at the free end of the beam to measure the response of the 
structure. An impact hammer was used to excite the beam with hammer hits made at various points 
along the length the beam. The first nine natural frequencies were estimated from the experimental 
data using a Nyquist plot circle fitting technique. These same natural frequencies were calculated 
by numerical solution of Eq. (28) for the first nine modes. These results are given in table four. 
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Apparently, the estimated frequencies agree very well with the theoretically predicted ones. In fact, 
the standard deviation of the estimated frequencies to the theoretical is generally small. The error 
between the estimates and predictions is on the order of 1%. 

SUMMARY 

Actuators that can be modeled as lumped second order systems were examined for use in vibration 
control of a distributed parameter system. A finite dimensional model provided insight that was 
extended to the distributed parameter case. It was seen that ignoring the actuator dynamics can lead 
to an unstable control law formulation. Secondly, the feedback paths available for output feedback 
control were identified and examined in terms of closed loop stability. This resulted in closed loop 
stability conditions for computing stable control gains. Finally, an example was given where it 
was seen that the use of noncolocated feedback gave better performance than solely colocated 
feedback. 

An infinite dimensional formulation of a cantilevered beam with actuator dynamics was presented 
and experimentally verified. This model included both structural damping, in the form of viscous 
and Kelvin-Voight damping, and the actuator dynamics. It remains to complete the computational 
studies of the infinite dimensional model and its approximations. Initial experimental verification 
showed good agreement between theoretically predicted and experimentally estimated natural 
frequencies. 
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APPENDIX: ACRONYMS 

DPS - Distributed Parameter System 

LVDT - Linear V ariab le Differential Transformer 

MDOF - Mul tiple Degree of Freedom 

MPMA - Multiple Proof Mass Actuator 

PMA - Proof Mass Actuator 

SDOF - Single Degree of Freedom 
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Table 1. Beam parameters 


Property 

Units 

Symbol 

Value 

Length 

m 

1 

1.0 

Moment of Inertia 

m4 

I 

1.64 x 10-9 

Area, cross section 

m 2 

A 

0.597 x 10- 3 

Young’s Modulus 

GPa 

E 

26.8 

Viscous damping 

N s nr 1 

ci 

1.75 

Kelvin-Voight damping 

kN s nr 1 

c 2 

20.5 

linear Density 

kg nr 1 

P 

1.02 


Table 2. Actuator mass parameters 


Property 

Unit 

Symbol 

Value 

Proof mass 

kg 

nip 

0.225 

Deadmass 

is 



0.730 


Table 3. Example gain settings 


Gain 

System 1 

System 2 

Ka (N nr 1 ) 

23.9 

211 

g (N s nr 1 ) 

0 

77.5 

ca + eg (Ns m* 1 ) 

1.0 

17.5 

Cost J 

0.183 

0.0710 
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Table 4. Measured and predicted frequencies 


Mode 

Analytical Freq. 
(Hz) 

Experimental Freq. 
(Hz) 

Std dev. 

1 


3.48 

0.127 


17.7 

17.8 

0.233 


22.4 

22.6 

0.0479 


29.3 

28.5 

0.916 


68.3 

67.8 

0.517 


133 

132 

0.995 


219 

217 

1.30 

8 

327 

324 

3.24 

9 

457 

452 

4.40 



Figure 1. Single degree of freedom stucture with proof mass actuator (SDOF/PMA) 
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Initial condition response systems: (1) and (2) 



System 1 
System 2 


Time (sec) 

Figure 2. Structural position response: system (1) vibration absorber 

system (2) active control with velocity feedback 

Initial condition response systems: (1) and (2) 


! 

* r-J 
t/5 

a 

3 



System 1 
System 2 


Time (sec) 

Figure 3. Actuator mass response: system (1) vibration absorber 

system (2) active control with velocity feedback 
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actuator i 


actuator 2 


Figure 4. Multiple degree of freedom structure with multiple actuators (MDOF/MPMA) 
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